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ABSTRACT 


In  the  frequency  domain  all  materials  can  be  described  electrically  by  their 
complex  pennittivity  ( a  )  and  penneability  ( // ).  These  constitutive  parameters  determine 
the  response  of  the  material  to  electromagnetic  (EM)  radiation.  The  precise  knowledge  of 
complex  permittivity  and  penneability  is  required  not  only  for  scientific  but  also  for 
industrial  applications.  Due  to  the  uncertainties  in  manufacturing  processes,  often  the  only 
way  to  find  a  material’s  parameters  is  to  measure  them. 

The  concept  of  metamaterials,  exhibiting  negative  permittivity  and  penneability,  is 
attracting  a  lot  of  attention.  Such  materials  are  also  termed  left-handed  materials  (LHMs). 

This  thesis  examines  several  methods  to  determine  the  effective  pennittivity  and 
permeability  of  both  normal  materials  and  metamaterials.  CST  Microwave  Studio  is  used 
to  model  the  materials  in  both  free  space  and  rectangular  waveguide  environments  to 
calculate  the  S-parameters  (Sn  and  S21)  from  which  the  constitutive  parameters  can  be 
extracted. 
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I.  INTRODUCTION 


A.  BACKGROUND 

In  the  frequency  domain  all  materials  can  be  described  electrically  by  their 
complex  permittivity  ( s )  and  permeability  ( // )  which  determine  their  response  to 
electromagnetic  (EM)  radiation.  The  precise  knowledge  of  these  constitutive  parameters 
is  required  not  only  for  scientific  but  also  for  industrial  applications.  Due  to  the 
uncertainties  in  manufacturing  processes,  often  the  only  way  to  find  a  material’s 
parameters  is  to  measure  them. 

In  theory,  from  reflection  and  transmission  measurements,  the  complex  //  and  s 
can  be  extracted.  However,  there  are  several  practical  problems  that  arise,  including 
ambiguities  at  thicknesses  that  are  multiples  of  XI 4.  There  are  even  more  extraction 
difficulties  when  the  material  under  test  is  a  metamaterial,  which  has  complicated 
frequency  dependent  material  properties. 

The  transmission  and  reflection  measurements  can  be  done  in  free  space  with  a 
large  sample  of  material  if  care  is  taken.  However,  using  large  samples  is  undesirable  when 
many  different  material  parameters  are  being  investigated.  It  is  more  efficient  to  take  the 
measurements  in  a  waveguide  using  small  samples.  This  is  referred  to  as  a  waveguide 
simulator,  and  it  has  been  used  for  decades  in  the  design  of  phased  array  elements  [1]. 

For  metamaterials,  a  retrieval  method  based  on  simulation  is  very  attractive.  The 
parameters  of  the  physical  materials  used  to  construct  the  metamaterials  are  known  (e.g., 
the  dielectric  substrates  and  conductors).  But  the  effective  parameters  are  what  is  needed 
for  designing  devices  with  metamaterials.  The  effective  parameters  are  the  average  or 
macroscopic  values  that  would  be  presented  to  an  incident  field.  If  the  effective  values  can 
be  found  from  simulation,  it  is  much  more  efficient  than  building  samples  for  testing. 

A  method  that  uses  reflection  measurement  has  been  described  by  McGinnis  [2]. 
For  most  cases  this  general  method  should  work  equally  as  well  as  the  methods  shown  by 
Sun  [3]  and  by  McGinnis  [4].  The  latter  method  is  more  appropriate  if  the  material  has  very 
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high  permittivity,  permeability  or  attenuation,  or  the  measurement  varies  from  sample  to 
sample.  A  newer  method  [5]  using  the  same  concept  as  [4]  but  with  slightly  different 
formulas,  was  used  in  this  thesis. 

Metamaterials  have  complex  permittivity  and  permeability  that  are  both  negative 
(double  negative,  DNG)  over  some  frequency  band  resulting  in  a  negative  index  of 
refraction.  This  unusual  property  of  metamaterials  causes  electromagnetic  waves  to  be 
refracted  at  a  negative  angle  relative  to  conventional  materials.  Various  methods  to  retrieve 
their  constitutive  parameters  using  complex  valued  reflection  (Sn)  and  transmission  (S21) 
coefficients  have  been  developed.  One  of  them  is  to  parametrically  fit  the  S-parameter 
measurement  data  to  known  pennittivity  (e.g.,  a  Drude  model)  and  penneability  (e.g.,  a 
Lorentz  model).  However,  this  method  makes  distinct  assumptions  about  the  form  of  the 
pennittivity  and  penneability.  Therefore,  in  this  thesis,  methods  that  directly  retrieve  the 
pennittivity  and  permeability  are  examined  and  programmed. 


B.  OBJECTIVE 

The  objective  of  this  thesis  is  to  develop  and  verily  methods  to  retrieve  the 
constitutive  parameters  of  planar  materials  in  free  space  and  waveguide.  A  wave  matrix 
formulation  [6]  and  formulas  from  [4]  are  developed  to  calculate  the  reflection  and 
transmission  coefficients  for  arbitrary  thicknesses  of  materials.  Free  space  and  waveguide 
environments  are  used  to  extract  the  pennittivity  and  permeability  from  the  S-parameters. 
Numerical  simulation  software,  Microwave  Studio  (MWS)  by  CST,  is  used  to  model  the 
materials  and  calculate  the  S-parameters  for  comparison  purposes.  The  extraction 
calculation  is  done  in  a  MATLAB  program.  To  validate  the  program,  many  test  cases  are 
presented  and  simulation  results  are  compared  to  published  data. 


C.  ORGANIZATION  OF  THESIS 

In  Chapter  II,  metamaterials  and  their  unusual  characteristics  are  discussed.  The 

potential  application  of  metamaterials  in  waveguide  and  other  applications  are  presented. 
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In  Chapter  III,  the  free  space  environment  retrieval  method  is  proposed  and 

discussed.  The  simulation  set  up  using  MWS  for  the  free  space  environment  is  presented. 
The  results  extracted  from  the  S-parameters  are  compared  with  published  data. 

In  Chapter  IV,  the  waveguide  environment  simulation  set  up  and  retrieval  methods 
are  presented.  The  extracted  permittivity  and  permeability  are  described  and  compared. 

Finally  in  Chapter  V,  the  results,  conclusions  and  suggestions  for  future  studies  are 
discussed. 
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II.  METAMATERIALS 


A.  INTRODUCTION 

Metamaterials,  also  called  left-handed  materials  (LHMs),  are  a  class  of  composite 
materials  artificially  built  to  exhibit  unusual  properties  not  found  in  nature.  The  history  of 
LHMs  can  be  tracked  back  to  Veselago’s  theoretical  hypothesis  in  1968  [7],  in  which  he 
demonstrated  that  the  LHMs  would  result  in  unusual  optical  phenomena  when  light  passed 
through  them.  In  2000,  the  UCSD  (University  of  California,  San  Diego)  group,  which  was 
following  the  work  by  Pendry  et  al.  [8-10],  demonstrated  the  first  left-handed  material  [11, 
12].  This  LHM  made  use  of  an  array  of  conducting,  nonmagnetic  elements  and  an  array  of 
conducting  continuous  wires  to  achieve  a  negative  effective  permittivity  and  negative 
effective  permeability  respectively. 


B. 

NEGATIVE  PERMITTIVITY  AND  PERMEABILITY 

In  general,  for  a  passive  medium,  assuming  a  eJC0t  time  dependence, 

and 

s  =  srs0=(sr'-j£r")s0 

(2.1) 

(2.2) 

where  s0  =8.85x10  12F/m  and  //0  =  4;rx  10  7H/m.  For  most  materials  //  and  s  are 
scalars,  however,  for  anisotropic  materials,  matrices  must  be  used: 


and 


where,  for  example 


D  =  [s]E 


b  =  [m]h 


f 


zy 


(2.3) 

(2.4) 


(2.5) 
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The  elements  stj  (i,j  =  x,y  or  z  )  are  the  complex  relative  values  that  completely 

describe  the  permittivity  of  the  material.  A  similar  matrix  can  be  defined  for  the 
permeability: 


[/A  =  Ao 


Myx 

,  Mzx 


Mxy 

\ 

Mxz 

Vyy 

Myz 

Mzy 

Azzy 

(2.6) 


Metamaterials  are  artificial  materials  that  are  designed  to  have  unique 
electromagnetic  properties  such  as  a  negative  index  of  refraction.  For  a  DNG  metamaterial, 
one  should  write  the  permittivity  and  penneability  as  [13]: 

(2-7) 


and 


(2.8) 


The  wavenumber  and  the  wave  impedance  can  be  described  as  follows: 

and 


/T 

£r 


s[Jk  VI 

VG  v 


(2.9) 


(2.10) 


The  index  of  refraction  for  a  DNG  material  can  be  obtained  from  the  formula: 


n  = 


VgV/t  =-Vkl|/6 


(2.11) 


The  different  combinations  of  real  permeability  and  permittivity  with  various  signs 
are  listed  in  Table  1. 
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Table  1. 

Sign  of  Index  of  Rel 

fraction 

For  a  DNG  material  the  negative  sign  flips  the  angle  of  refraction  in  Snell’s  law. 
For  example,  a  negative  angle  of  refraction  equal  and  opposite  to  the  standard  angle  of 
incidence  can  be  clearly  observed  in  Figure  1. 


Figure  1 .  Electric  Field  Intensity  Distribution  for  the  Wave  Interaction  with  a  DNG 

Metamaterial  (After  Ref.  [13]) 

As  can  be  seen  in  Figure  2,  on  the  left,  a  ray  enters  a  negatively  refracting  medium 
and  is  bent  the  opposite  way  relative  to  the  surface  normal  compared  to  a  conventional 
right  handed  material.  On  the  right  is  the  wave  vector.  Negative  refraction  requires  that  the 
wave  vector  (direction  of  phase  fronts)  and  Poynting  vector  (direction  of  energy  flow) 
point  in  opposite  directions  [13]. 
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Figure  2.  Depiction  of  Negative  Refraction  (After  Ref.  [13]) 

The  traditional  way  to  build  a  metamaterial  with  both  negative  penneability  and 
permittivity  is  to  use  a  split  ring  resonator  (SRR)  and  wire  in  a  unit  cell.  The  metamaterial 
model  built  by  D.  R.  Smith  is  shown  in  Figure  3.  Macroscopically  this  material  behaves  as 
a  solid  material  with  an  effective  permittivity  and  penneability.  This  is  not  surprising  if  the 
cell  is  much  smaller  than  the  wavelength  of  the  wave  interacting  with  the  material. 
Flowever  it  has  been  found  that  cell  sizes  on  the  order  of  a  fifth  of  a  wavelength  are  small 
enough  for  the  material  to  behave  as  a  homogeneous  slab. 


Figure  3.  A  Negative  Index  Metamaterial  Formed  by  SRRs  and  Wires  Deposited  on 
Opposite  Sides  of  a  Standard  Circuit  Board.  (From  Ref.  [14]) 
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There  are  also  new  concepts  for  DNG  metamaterials.  In  Chapter  III  and  Chapter 
IV,  these  “new”  metamaterials  were  investigated,  one  of  which  uses  short  wire  pairs  while 
the  other  uses  H-type  pairs  as  inclusions.  The  major  difference  between  “traditional”  and 
“new”  metamaterials  is  that  “new”  metamaterials  are  relatively  easy  to  fabricate,  which 
make  them  easier  to  use  and  may  hasten  their  development. 


C.  APPLICATIONS  OF  METAMATERIALS 

Metamaterials  have  attracted  a  lot  of  attention  in  recent  years.  Some  of  the  possible 
applications  of  metamaterials  are  discussed  below. 

1.  Waveguide  Miniaturization 

The  physical  size  of  a  waveguide  is  an  important  practical  parameter  and  there  is 
an  increased  need  for  miniaturized,  compact  waveguides.  In  [15],  a  method  to  minimize 
the  waveguide  using  uniaxial  negative  penneability  metamaterial  was  proposed.  By 
putting  SRRs  into  the  waveguide,  a  backward  passband  can  occur  below  the  cut-off 
frequency  of  the  waveguide.  In  this  case,  the  waveguide  supports  propagation  of 
backward  waves  in  this  passband  below  the  cut-off  frequency  of  the  dominant  TE  mode. 
The  characteristic  for  this  combination  presents  a  chance  to  construct  a  smaller 
waveguide,  but  still  allow  the  electromagnetic  waves  to  propagate.  The  schematic  of  the 
waveguide  loaded  with  SRRs  is  as  shown  in  Figure  4. 
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Figure  4.  (a)  The  Waveguide  Filled  with  the  Split  Ring  Resonators  (b)  The  Unit  Cell 

(From  Ref.  [15]) 
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The  difference  between  the  waveguides  loaded  with  SRRs  with  the  resonant 
frequency  of  SRRs  above  the  cut-off  frequency  fc  and  below  /’  (  /’  =  1 2.5  GHz  for 
<7  =  12  mm,  fc  =4.28  GHz  for  a  =  35  mm)  is  described  in  Figure  5.  Note  that  when 
the  resonant  frequency  of  SRRs  (/0)  is  higher  than  the  cut-off  frequency  of  the 

waveguide,  a  stop  band  with  negative  transversal  permeability  is  observed.  On  the  other 
hand,  when  f0<  fc,  a  backward  passband  is  observed  below  the  cut-off  frequency  of  the 
waveguide. 
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Figure  5.  (a)  Calculated  Power  Transmission  Coefficient  of  the  Waveguide  (a= 35  mm) 

Filled  with  Metamaterial  Based  on  Resonant  Inclusions 
( f0  =  7.8  GHz ,  fmp  =  7.9  GHz ,  y  =  1  MHz  ,  and  f0  >  fc ).  (b)  Calculated  Power 
Transmission  Coefficient  of  the  Waveguide  (a=  12  mm)  Filled  with  Metamaterial 
Based  on  Resonant  Inclusions  (/0  =  7.8  GHz ,  fmp  =  7.9  GHz ,  y  =  1  MHz  ,f0<  fc ). 

(From  Ref.  [15]) 
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2.  Directional  Coupler 

A  DNG-based  guiding  structure  offers  a  possibility  for  the  “backward”  coupling 
between  such  guides  and  the  conventional  ones.  If  one  of  the  coupled  waveguides  carries 
power  in  one  direction,  the  second  waveguide,  placed  in  its  vicinity,  will  send  some  of 
this  power  in  the  opposite  direction.  As  the  distance  between  the  two  guides  decreases, 
the  coupling  efficiency  of  the  redirected  power  increases  exponentially.  A  schematic  of 
this  kind  of  coupler  is  illustrated  in  Figure  6. 
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Coupled  Out 
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Figure  6.  Schematic  of  a  Metamaterial-based  Directional  Coupler  in  Microstrip 

Implementation  (From  Ref.  [16]) 


This  kind  of  DNG-based  directional  coupler  has  been  shown  to  have  improved 
coupling,  a  wider  range  of  coupling  values  than  in  conventional  implementations,  and  a 
smaller  coupling  length. 

An  experiment  has  been  performed  by  Eleftheriades  [17]  where  he  uses  two  sets 
of  couplers  to  compare  the  coupled  power  and  isolation  between  Microstrip 
/Negative-Refractive-Index  (MS/NRI)  and  Microstrip/Micro  strip  (MS/MS).  The  results 
are  shown  in  Figures  7  through  10. 
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Figure  7.  Microstrip/Negative-Refractive-Index  and  Microstrip/Microstrip  Couplers  of 
Equal  Length,  Line  Spacing,  and  Propagation  Constants  (Lrom  Ref.  [17]) 
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Ligure  8.  A  Comparison  of  the  Coupled  Power  Levels  for  the 
Microstrip/Negative-Refractive-Index  and  the  Microstrip/Microstrip  Couplers 

(Lrom  Ref.  [17]) 
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Si 4  ;  solid  line  -  experiment ;  broken  line  -  simulation 


Figure  9.  A  Comparison  of  the  Isolation  for  the  Microstrip/Negative-Refractive -Index 

and  the  Microstrip/Microstrip  Couplers  (From  Ref.  [17]) 


Figure  10.  The  Physical  Mechanism  by  which  Power  Continuously  Leaks  from  the 
Microstrip  Line  to  the  Negative -Refractive-Index  Line  (From  Ref.  [17]) 
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For  MS/NRI,  the  coupled  power  is  6  dB  higher  than  that  in  MS/MS  (Figure  7)  and 
the  isolation  for  MS/NRI  is  about  4  dB  lower  than  that  in  the  MS/MS  structure  (Figure  8). 
The  new  MS/NRI  coupler  exhibits  better  performance  than  the  traditional  MS/MS 
coupler  in  terms  of  higher  coupling  power,  isolation,  and  lower  return  loss.  The  coupling 
mechanism  is  depicted  in  Figure  10. 

There  are  more  applications  that  can  exploit  the  unique  characteristics  of 
metamaterials,  such  as  [16]: 

1 .  Wideband  directive  antennas 

2.  Perfect  lens  and  superlens 

3.  Phase  compensators  and  sub-wavelength  resonators 

4.  Guided  wave  applications 

5.  Electrically  small  antennas 

6.  Filtering  and  beam  shaping 

The  two  examples  and  the  list  above  include  only  a  few  from  a  wide  variety  of 
applications.  As  more  and  more  people  become  aware  of  these  new  materials  and  their 
properties,  one  can  expect  to  see  more  applications  using  metamaterials.  In  the  next 
chapter  the  extraction  of  //  and  s  for  bulk  materials  in  free  space  is  discussed.  The 
procedure  is  applied  to  both  conventional  materials  and  metamaterials. 
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III.  FREE  SPACE  ENVIRONMENT 


A.  INTRODUCTION 

In  this  chapter,  measurements  on  materials  in  the  free  space  environment  are  used 
to  extract  the  constitutive  parameters  from  scattering  parameters.  In  order  to  verify  the 
accuracy  of  scattering  parameters  generated  from  simulation  using  MWS,  the  wave  matrix 
approach  [18]  is  used  to  calculate  the  reflection  and  transmission  coefficients  and  compare 
with  those  generated  from  simulation. 


B.  FORMULATION 

The  uniform  plane  waves  traveling  in  the  z  direction  through  an  infinite  slab  with 
relative  permittivity  sr  ,  penneability  //.  ,  wave  impedance  //,  and  thickness  t  is 
considered.  Figure  1 1  illustrates  the  simplified  model  described  above. 
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Figure  11.  Schematic  of  Normal  Incidence  at  an  Infinite  Slab  in  Free  Space 

The  general  form  of  the  wave  matrix  formula  for  a  three-layer  medium  that  was 
used  to  examine  the  accuracy  of  the  simulation  results  is  [18] 
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If  the  thickness  of  the  third  layer  extends  to  go  ,  then  it  is  possible  to  state  that  b3=  0  .  The 
total  transmission  and  reflection  coefficients  of  the  slab  become 


T  = 


C3 
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C1 


(3.2) 


and 


r  _  b\_  _  a2\ 


(3.3) 


The  difference  between  the  calculated  and  simulated  transmission  and  reflection 
coefficients  for  a  normal  material  (both  selected  //  and  s  are  positive)  are  presented  in 
the  next  section. 

The  retrieval  method  used  in  the  first  part  of  this  thesis  is  based  on  the  approach 
used  by  Xudong  Chen  et  al  in  [19],  where  the  S-parameters  are  defined  in  terms  of  the 
reflection  and  transmission  coefficients  as 


Sn=R 


(3-4) 


and 


S2l=Te 


jk0d 


(3.5) 


where  k0  is  the  wave  number  in  free  space  and  d  is  the  thickness  of  the  material  being 
examined.  The  impedance  Z  and  the  refractive  index  n  are 


Z  =  ± 


(\  +  Su)2-S2; 


]l(i-sny-s212 


(3.6) 


and 


n  = 


k0d 


Im 


ln(e"'v)' 


+  2ni7t 


Re 


ln(e"’v) 


(3.7) 


18 


where  m  is  an  integer  related  to  the  branch  index  of  the  real  part  of  n  and  the  value  of 
ejnk"d  is  obtained  from  Su,  S2] ,  and  z  from 


]nk0d  _ 
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z+T 


(3.8) 


As  described  in  [19],  Z  and  n  are  related  so  that  their  relationship  can  be  used  to 
determine  the  correct  sign  of  Z  .  Two  cases  are  distinguished  in  order  to  correctly  find  the 
sign  of  Z  .  The  first  case  is  for  |Re(Z)|  >  8  ,  where  8  is  a  positive  number  and  for  which 
Re(Z)  >  0  .  In  the  second  case,  the  sign  of  Z  is  chosen  so  that  the  corresponding 


\m(n)  >  0 ,  or  equivalently 
Re(«)  can  be  found  in  [19]. 


,y«v 


<  1 .  The  method  to  precisely  determine  the  branch  of 


C.  RETRIEVAL  USING  NUMERICALLY  GENERATED  S-PARAMETER 

In  the  first  part  of  this  section,  it  is  demonstrated  that  the  simulated  data  is  accurate 
by  comparing  it  with  the  data  calculated  from  the  wave  matrix  formula.  The  effectiveness 
of  the  retrieval  program  is  also  examined  by  extracting  the  constitutive  parameters  of  the 
model  and  comparing  them  with  the  original  values  used  to  compute  T  and  R,  and 
published  paper  data. 

Simulations  were  done  with  CST  Microwave  Studio  (Computer  Simulation 
Technology  GmbH,  Darmstadt,  Germany),  which  uses  the  finite-integration  method  in  the 
time  domain  to  detennine  reflection  and/or  transmission  properties  of  any  given  structures. 

1.  Model  1  -  Normal  Material 

The  model  used  here  is  a  cube  with  dimension  2.5  mm.  The  material  is  given 
the  constitutive  parameters  =3  and  sr  =5.  Boundary  conditions  are  set  up  to  be 

magnetic  (Ht=  0)  in  the x  direction,  electric  (Et  =0)  in  the  v  direction,  and  open  in 

the  z  direction;  therefore  an  infinite  slab  can  be  modeled  as  shown  in  Figure  12.  The 
plane  wave  propagates  in  the  negative  z  direction. 
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Figure  12.  Normal  Material  with  nr  =  3  and  £,.  =  5 ,  and  2.5  mm  side  length 


The  differences  between  the  calculated  (using  wave  matrix  fonnula)  and 
simulated  (using  MWS)  transmission  ( T )  and  reflection  (  R  )  coefficients  are  shown  in 
Figures  13  and  14. 


difference  between  data  and  wave  matrix  difference  between  data  and  wave  matrix 


Frequency,  GHz 


Frequency,  GHz 


Figure  13.  (a)  Magnitude  Differences  of  Transmission  Coefficient  between  the 

Calculated  and  Simulated  Data  (b)  Phase  Differences  (Degrees)  of  Transmission 
Coefficient  between  the  Calculated  and  Simulated  Data 
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difference  between  data  and  wave  matrix  difference  between  data  and  wave  matrix 


Frequency,  GHz  Frequency,  GHz 


Figure  14.  (a)  Magnitude  Differences  of  Reflection  Coefficient  between  the  Calculated 

and  Simulated  Data  (b)  Phase  Differences  (Degrees)  of  Reflection  Coefficient 
between  the  Calculated  and  Simulated  Data 

According  to  Figures  13  and  14,  the  differences  between  the  calculated  and 
simulated  data  are  almost  zero;  therefore  it  can  be  concluded  that  the  software  set  up  is 
correct  and  an  infinite  panel  is  simulated.  After  this  simulation  accuracy  verification,  the 
constitutive  parameter  retrieval  procedure  can  be  conducted  based  on  the  same  set  up. 

For  the  same  model  and  data,  an  extraction  program  was  used  to  retrieve  the 
permittivity  and  penneability  from  the  scattering  parameters  generated  by  MWS.  The 
results  of  the  extracted  //,.  and  sr  using  MATLAB  are  shown  in  Figures  15  and  16. 
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Figure  15. 


Figure  16. 
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Extracted  Real  Part  and  Imaginary  Part  of  Permeability  of  Model  1 
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Extracted  Real  Part  and  Imaginary  Part  of  Permittivity  of  Model  1 
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The  retrieval  program  correctly  extracts  the  relative  permeability  (5)  and 
permittivity  (3)  from  the  simulation  data.  Thus,  the  retrieval  program  works  fine  with 
nonnal  RH  material  whose  constitutive  parameters  are  both  positive.  Next,  the  process  is 
extended  to  LH  materials. 

2.  Model  2  -  “Traditional”  Metamaterial 

The  “traditional”  metamaterial  model,  consists  of  a  split  ring  resonator  (SRR)  and 
wire  together  to  form  a  material  with  both  negative  permittivity  and  permeability  over 
some  frequencies  (see  Figure  17).  The  dimensions  of  this  model  are  the  same  as 
described  in  [20].  The  unit  cell  is  cubic,  with  a  cell  dimension  of  d  =  2.5  mm.  A 
0.25  mm  thick  substrate  with  s  =  4.4,  loss  tangent  of  0.02  is  assumed.  The  copper 
thickness  is  0.017  mm.  The  width  of  the  wire  is  0.14  mm,  and  it  runs  the  length  of  the 
unit  cell.  The  outer  ring  length  of  the  SRR  is  2.2  mm  and  both  rings  have  a  linewidth  of 
0.2  mm.  The  gap  in  each  ring  is  0.3  mm,  and  the  gap  between  the  inner  and  outer 
rings  is  0.15  mm.  The  schematic  of  this  metamaterial  slab  is  shown  as  Figure  18. 
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Figure  17.  Model  2  Built  in  MWS  with  SRR  and  Wire.  The  SRR  Provides  Negative 
Permeability  While  the  Wire  Provides  Negative  Permittivity  Just  Above  the 
Resonant  Frequency  of  the  Material 


23 


Figure  18. 


Schematic  of  the  Metamaterial  Slab  with  Uniform  Plane  Wave  Traveling  in 
the  Negative  z  Direction  (Five  Layers  of  Unit  Cells  in  Both  x  and  y 

Directions  are  Shown  Here) 


The  same  set  up  as  for  Model  1  was  used,  which  has  the  boundary  condition 
Ht-  0  in  the  x  direction,  Et  =  0  in  the  y  direction,  and  open  in  the  z  direction.  The 

boundary  condition  of  this  kind  not  only  ensures  that  the  magnetic  field  goes  through 
SRR  but  also  ensures  that  the  model  represents  an  infinite  slab. 

The  extracted  permittivity,  permeability,  index  of  refraction,  and  impedance  are 
shown  in  Figures  19  through  22. 
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Figure  19.  Extracted  Permittivity  of  Model  2 


Figure  20.  Extracted  Permeability  of  Model  2 
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Figure  2 1 .  Index  of  Refraction  of  Model  2 


Figure  22.  Impedance  of  Model  2 
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The  extracted  constitutive  parameters,  index  of  refraction,  and  impedance  of  the 
material  are  very  close  to  the  published  data  [20].  However,  in  the  simulation  process,  it 
was  found  that  the  meshing  property  is  very  important  in  terms  of  the  accuracy  of  the 
scattering  parameters  and  the  position  of  the  resonant  frequency.  Because  the  author  does 
not  have  the  exact  meshing  parameters  which  the  published  paper  used,  the  final  results 
are  not  identical  to  the  published  data,  although  they  are  very  close. 

3.  Model  3  -  “New”  Metamaterial  Consisting  of  Short  Wire  Pairs 

The  third  model  that  was  used  to  verify  the  effectiveness  of  the  program  is  based 
on  the  same  structure  proposed  in  [21].  The  dielectric  layer  has  a  thickness  of  254  //m 
and  a  relative  dielectric  constant  of  2.53.  The  copper  wire  pairs  have  a  thickness  10  //m. 
The  widths  of  all  metal  lines  are  1  mm.  The  length  of  the  short  wire  pairs  was  7  mm 
and  the  unit  cell  size  was  9.5x7x2.274  mm3.  In  the  structure  arrangement,  the 
conventional  SRR  is  replaced  with  a  pair  of  short  parallel  wires  and  the  continuous  wires 
are  preserved  (see  Figure  23).  This  “new”  structure  of  metamaterial  makes  the  fabrication 
procedure  much  easier.  The  short  wire  pairs  give  the  negative  permeability  and  the 
continuous  wires  give  the  negative  permittivity  over  certain  frequencies.  The  schematic 
of  this  metamaterial  slab  is  shown  as  Figure  24. 


Figure  23. 
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Figure  24.  Schematic  of  the  Metamaterial  Slab  with  Uniform  Plane  Wave  Traveling  in 
the  Negative  z  Direction  (Six  Layers  of  Unit  Cells  in  Both  x  and  y 
Directions  are  Shown  Here) 


Figures  25  through  27  show  the  extracted  permeability,  permittivity  and  the 
refractive  index  using  the  simulation  data  and  program.  Note  that  the  simulation  result  is 
very  sensitive  to  the  setting  of  “mesh  cells”  in  MWS.  Different  mesh  cells  can  generate 
slightly  different  resonant  frequencies.  One  way  to  avoid  this  uncertainty  is  to  use  the 
function  “adaptive  mesh  refinement”  in  MWS.  This  function  is  able  to  calculate  the 
S-parameters  with  different  mesh  cells  until  the  difference  between  the  S-parameters  is 
less  than  the  specified  value.  This  method  can  ensure  that  the  mesh  size  is  small  enough 
to  make  the  results  stable. 
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Figure  25.  Extracted  Permeability  of  Model  3 


Figure  26.  Extracted  Permittivity  of  Model  3 
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Figure  27.  Refractive  Index  of  Model  3 

In  comparing  the  figures  with  published  data  in  [21],  it  can  be  seen  that  the 
program  discussed  in  this  thesis  for  the  constitutive  parameter  retrieval  from  scattering 
parameters  not  only  can  be  used  for  “traditional”  metamaterials,  but  also  for  the  “new” 
metamaterials. 

4.  Model  4  -  “New  Metamaterial”  Consisting  of  H  Type  Inclusions 

A  fourth  model  considered  is  also  a  new  metamaterial  with  H  type  inclusions. 
This  model  is  built  based  on  the  same  dimensions  described  in  [22].  One  unit  cell  of  the 
H  type  structure  is  shown  in  Figure  28.  The  dielectric-layer  thickness  is  254  //m  and  the 
dielectric  constant  is  2.53.  The  copper  was  formed  in  the  wire  pairs  patterns  have  10  jura 
thickness.  The  width  of  the  metal  wires  was  1  mm.  The  length  and  width  of  the  metal 
bars  at  each  end  of  the  wire  pairs  were  4  and  0.2  mm,  respectively.  The  length  of  the 

short  wire  pairs  was  4  mm,  and  the  unit  cell  size  was  4x8x2.274  mm3.  The 
schematic  of  this  metamaterial  slab  is  shown  as  Figure  29. 
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z 

Figure  28.  One  Unit  Cell  of  Model  4  Consists  of  Two  H  Type  Inclusions 


x 


Figure  29.  Schematic  of  the  Metamaterial  Slab  with  Unifonn  Plane  Wave  Traveling  in 
the  Negative  z  Direction  (Six  Layers  of  Unit  Cells  in  Both  x  and  y 
Directions  are  Shown  Here) 


In  order  to  test  the  sensitivity  of  this  kind  of  material  relative  to  the  dimension  of 
the  inclusions  and  the  number  of  mesh  cells  being  used  for  the  simulation,  the  following 
two  scenarios  were  examined. 

a.  Change  the  Width  of  the  Bars  at  Each  End  of  the  Wire  Pairs 
The  width  of  the  metal  bars  at  each  end  of  the  wire  pairs  was  slightly 
changed  in  order  to  see  the  influence  on  the  S-parameters.  The  width  was  changed  from 
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0.2  min  to  0.25  mm  and  the  adaptive  mesh  refinement  function  was  enabled  so  as  to 
obtain  optimal  results.  The  S-parameters  for  the  original  0.2  nun  and  the  modified  0.25 
mm  structures  are  presented  in  Figures  30  and  31. 


|S2, 1 1  versus  Passes 


Figure  30.  Magnitude  of  S21  of  the  Original  0.2  mm  Model 


|S2, 1 1  versus  Passes 


Figure  3 1 .  Magnitude  of  S21  of  the  Modified  0.25  mm  Model 


The  resonant  frequency  in  Figure  30  with  the  optimized  S2l  value  occurs 
around  15.65  GHz  (purple  curve).  The  optimized  resonant  frequency  in  Figure  31  (width 
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changed  by  0.05  mm)  shifted  to  the  left  to  around  15.25  GHz  (purple  curve).  This 
observation  shows  that  the  resonant  frequency  is  very  sensitive  to  the  dimension  of  the 
inclusions. 

b.  Using  More  Mesh  Cells  at  Simulation  Starting  Point 
For  the  purpose  of  understanding  the  influence  of  using  a  different  number 
of  mesh  cells  at  the  simulation  starting  point,  the  meshing  parameters  in  MWS  were  set 
up  to  make  the  total  number  of  mesh  cells  higher  than  the  amount  used  in  scenario  1 
shown  in  Figure  30.  The  simulation  result  is  shown  in  Figure  32. 


Figure  32.  Magnitude  of  S21  of  the  Same  Model  in  Figure  29  Except  for  Using  More 

Mesh  Cells  as  the  Simulation  Starting  Point 


It  can  be  clearly  seen  that  using  more  mesh  cells  at  the  starting  point 
(Figure  32)  has  the  resonant  frequency  of  the  optimized  S2 ,  at  about  15.75  GHz  (blue 
curve).  On  the  other  hand,  the  one  using  less  mesh  cells  (Figure  30)  has  the  resonant 
frequency  of  the  optimized  S21  at  about  15.65  GHz.  This  is  due  to  the  increase  in 

discretization  and  simulation  errors  when  the  mesh  size  is  not  optimized  and  is  either  too 
large  or  too  small.  As  a  result,  the  starting  point  of  the  mesh  cells  number  should  also  be 
considered  in  order  to  get  the  real  optimal  value.  The  other  way  to  make  sure  the  real 
value  is  obtained  is  to  compare  with  measurements. 

To  compare  the  retrieved  results  with  the  published  paper,  the  closest  S-parameter 
offered  in  the  paper  [22],  which  is  shown  in  Figure  30  pass  2  (green  curve),  was  chosen. 
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The  extracted  permeability,  permittivity  and  the  index  of  refraction  are  provided  in  Figures 
33  through  36  and  the  results  from  [22]  are  in  Figure  37  and  38  for  comparison  purposes. 


|S2,1|  versus  Passes  ri  CB 


Frequency  /  G-fc 

Figure  33.  Simulated  S21  of  the  H  type  Inclusions  Metamaterial  That  is  Closest  to  the 

Published  Data 


Figure  34.  Extracted  Permeability  of  Model  4 
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F igure  3  5 .  Extracted  Permittivity  of  Model  4 


Figure  36.  Refractive  Index  of  Model  4 
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Frequency  (GHz) 


Figure  37.  (a)  Extracted  Permittivity  e  and  (b)  Permeability  // ,  Using  the  Simulated 

(Red  Solid  Curves)  and  Measured  (Blue  Dotted  Curves)  Transmission  and 
Reflection  Data  (From  Ref.  [22]) 


Figure  38.  Extracted  Refractive  Index  n  ,  Using  the  Simulated  (Solid  Curves)  and 

Measured  (Dotted  Curves)  Transmission  and  Reflection  Data.  The  Red  and  Blue 
Curves  Show  the  Real  Part  of  n  and  Imaginary  Part  of  n  ,  Respectively  (From 

Ref.  [22]) 
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The  extracted  constitutive  parameters  and  the  refractive  index  are  very  close  to  the 
published  data.  Again  the  discrepancy  is  caused  by  the  simulation  set  up  which  is  very 
sensitive  to  the  dimension  of  inclusions  and  the  number  of  mesh  cells  used. 

D.  SUMMARY 

In  this  chapter,  four  models  were  used  to  extract  the  permeability  and  pennittivity 
from  the  simulated  S-parameters.  The  formula  used  was  based  on  the  free  space 
environment;  therefore  boundary  conditions  had  to  be  defined  so  that  the  model  could  be 
treated  as  an  infinite  slab.  Test  cases  were  performed  to  verify  the  effectiveness  of  the 
simulation  set  up.  It  was  found  that  the  differences  between  the  simulated  and  calculated 
S-parameters  were  negligible.  The  resonant  frequencies  of  metamaterial  models  were  very 
sensitive  to  both  the  dimensions  and  mesh  cell  number;  therefore,  it  is  good  to  have  both 
simulation  and  measurement  data  for  comparison  purposes.  The  above  models  have  shown 
that  the  extraction  program  discussed  in  this  thesis  can  work  not  only  for  normal  materials 
but  also  metamaterials. 
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IV.  WAVEGUIDE  ENVIRONMENT 


A.  INTRODUCTION 

In  this  chapter,  a  retrieval  method  for  materials  in  a  waveguide  environment  is  used 
to  extract  the  constitutive  parameters  of  the  materials.  The  purpose  of  this  waveguide 
method  is  to  demonstrate  that  the  method  of  images  can  be  applied  to  structures  for 
permeability  and  pennittivity  retrieval.  Therefore  a  small  sample  of  material  is  enough  to 
perform  the  extraction  or  measurement  procedure.  The  results  from  the  free  space 
environment  method  are  also  compared  to  support  this  hypothesis. 

The  image  principle  says  that  if  an  object  in  a  metallic  waveguide  can  be  imaged 
properly  with  respect  to  the  walls,  after  removing  the  waveguide  walls,  the  object  can  be 
treated  as  an  infinite  structure  in  unbounded  space.  The  equivalent  structure  is  composed  of 
an  infinite  number  of  repeated  objects  in  the  directions  perpendicular  to  the  waveguide 
walls.  An  example  of  the  image  principle  applied  to  the  waveguide  model  is  shown  in 
Figure  39. 
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y 
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B.  FORMULATION 

If  the  imaging  is  done  properly,  waves  propagating  down  the  waveguide  appear  to 
be  incident  onto  an  infinite  panel.  Generally  only  the  dominant  TEW  mode  is  allowed  to 
propagate.  It  can  be  decomposed  into  two  plane  wave  components 


f  7TX  ^ 

j[  «■--/?*)  -/ 

E_  =  E0  cos 

=  2  E 

eV  a  )  +  e  V  a  ) 

l  a  J 

(4.1) 


and  the  angles  of  incidence  with  respect  to  the  sample  face  normal  are  ±  tan 

Only  a  few  angles  will  be  possible  given  the  restrictions  on  imaging  the  structure.  Note  that 
the  electrical  size  of  a  must  be  large  to  obtain  near  normal  incidence  but  should  be  limited 
to  only  have  TEl0  mode.  More  cells  must  be  used,  which  increases  the  number  of  mesh 
nodes  and  simulation  time. 

Consider  a  material  of  thickness  8  that  is  placed  in  a  waveguide.  The  distances 
from  the  front  face  and  back  face  of  the  material  to  the  two  openings  of  the  waveguide  are 
d]  and  d2  .  The  material  is  perpendicular  to  the  direction  of  the  propagating 

electromagnetic  waves.  It  is  assumed  that  only  the  fundamental  waveguide  mode  TEi0 
can  propagate.  The  geometry  for  the  waveguide  retrieval  method  is  shown  below  in  Figure 
40.  The  scattering  parameters  can  be  determined  by  measurement  with  a  physical  system, 
or  by  simulation  with  computational  electromagnetic  (CEM)  software. 
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Figure  40.  Schematic  of  Waveguide  S-parameter  Retrieval  Setup  (After  Ref.  [4]) 


The  formulas  used  to  verify  the  MWS  model  set  up  are  from  [4].  In  the  case  of  a 
given  /ur  and  sr,  the  S-parameters  Sn(  =  Rm)  and  S21  ( =  Tm  )  can  be  obtained  from  the 
equations  that  follow. 


First,  the  propagation  constant  in  the  waveguide  filled  with  free  space  and  the 
material  are  defined  as  /?(l  and  J32  respectively, 


and 


2 

yC  ) 

(4.2) 


Pi  =  Mr£r 
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(4.3) 


The  reflection  coefficient  and  transmission  coefficient  at  the  material  surfaces  can 
be  described  as  follows: 


R  =  R2 


1  -e 


.Mid 


R2  -  e 


np2s 


(4.4) 
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and 


T  =  e~iPlS  RR2) 


(4.5) 


with 


=  mA>zA 

Mrfio  +  Pi 


(4.6) 


For  the  case  where  the  waveguide  is  not  fully  filled,  which  is  the  case  in  Figure  40, 
the  distances  between  material  surfaces  and  the  openings  of  the  waveguide  d]  and  d2 
must  be  included  in  the  formula  in  order  to  have  the  correct  phase  of  the  S-parameters. 
Therefore,  the  final  reflection  and  transmission  coefficients  can  be  described  as 

Sn=Rm=Re~j2M  (4.7) 


and 


S*=Tm=Te 


Jflo{d\+d  2) 


(4.8) 


This  formula  was  compared  to  the  MWS  simulated  S-parameters  to  make  sure  the 
simulation  set  up  in  the  waveguide  was  correct.  The  comparison  is  shown  in  the  next 
section. 

The  retrieval  method  used  in  the  second  part  of  this  thesis  is  based  on  the  approach 
used  by  Zwick  et  al  in  [23],  where  the  S-parameters  are  defined  in  terms  of  the  reflection 
and  transmission  coefficients  as 
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R(l~T2)  j2rA 
1  -r2t2 


s 
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r(i-i?2) 

1  -r2t2 


e~jro(d\+d2) 


Equations  (4.9)  and  (4.10)  can  be  solved  for  R  : 

R_K±^K2-4Sn2 
2  Sn 


(4.9) 

(4.10) 


(4.11) 
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with 


i:  =  sn2-s212+i 


(4.12) 


Sn+S2l-R 

1  -r{su+s21) 


(4.13) 


In  Equation  (4.1 1),  the  sign  of  the  square  root  has  to  be  chosen  such  that  \R\  <  1 . 
The  parameters  y  and  y0  can  be  calculated  as  follows: 


J2kft-ln(T) 


(4.14) 


7o  =  y'2;r,  f-ju0s0  -  — 
V  V  2# 


(4.15) 


with  k  being  an  integer  value,  /  the  frequency,  and  a  the  width  of  the  waveguide.  The 
details  on  how  to  choose  the  correct  k  is  described  in  [23],  Basically,  k  is  first  split  into 
a  frequency-dependent  part  k  (/)  and  an  offset  k0  as  follows: 


k  =  k0  +  k(f) 


(4.16) 


where  k(  f  )  is  simply  determined  by  using  the  fact  that  the  phase  of  T  has  to  be 
monotonically  decreasing  without  large  steps.  The  constant  k(]  can  be  obtained  by  finding 
the  k0  for  which  sr  and  //,  are  almost  constant  over  frequency.  The  complex 
permeability  and  pennittivity  are  finally  determined  as 


7(1  +  *) 
^07o(l-*) 


(4.17) 
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C.  RETRIEVAL  USING  NUMERICALLY  GENERATED  S-PARAMETERS 
1.  Model  5  -  Normal  Material  in  Waveguide 

Model  5  is  a  normal  (RH)  material  slab  with  relative  permeability  (3)  and 
permittivity  (5).  The  dimension  of  this  slab  is  8x4.2x2.274  mm3.  Boundary  conditions 
are  chosen  to  be  electric  (El  =0)  on  the  top  and  bottom  and  side  walls.  The  schematic 
of  this  model  is  shown  in  Figure  41. 


Figure  41.  Normal  Slab  Material  with  jur=  3  and  £,.=5  in  a  Waveguide.  The 

Dimension  of  the  Slab  is  8x  4.2  x  2.274  mm3. 

The  operating  frequency  must  be  above  the  cut-off  frequency  of  the  waveguide 
such  that  the  electromagnetic  waves  can  propagate  through  the  waveguide.  In  the 
simulation,  it  is  required  that  only  the  fundamental  TEl0  mode  can  propagate.  Therefore, 

not  only  is  it  necessary  to  operate  above  the  cut-off  frequency  but  it  is  also  necessary  to 
prevent  higher  modes  from  being  generated.  Cut-off  frequency  for  the  waveguide  can  be 
calculated  from  the  formula 
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For  this  model,  the  width  of  the  waveguide  is  a  =  8  nun  and  the  height  of  the 
waveguide  is  b  =  4.2  mm,  and  //  =  //„  and  £  =  s0  for  free  space.  The  operating 
frequency  range  for  this  waveguide  can  be  chosen  from  Table  2. 


(fc) 

\  ^  6  /  mn 

m  =1 

m  =2 

m  =3 

n  =0 

18.74  GHz 

37.48  GHz 

56.22  GHz 

n  =1 

40.32  GHz 

51.76  GHz 

66.60  GHz 

n  =2 

73.81  GHz 

80.64  GHz 

90.88  GHz 

n  =3 

108.72  GHz 

113.46  GHz 

120.96  GHz 

Table  2.  Cut-off  Frequencies  for  Waveguide  with  a  =  8  mm  and  b  =  4.2  mm  . 


According  to  Table  2,  the  operating  frequency  should  be  chosen  between  18.74 
and  37.48  GHz  to  only  allow  the  TEW  mode  to  propagate. 

As  described  in  Section  C  of  Chapter  III,  the  calculated  and  simulated 
transmission  (T)  and  reflection  (7?)  coefficients  were  compared  in  order  to  verify  the 

set  up  in  MWS.  The  formulas  used  were  the  ones  described  in  Section  B.  The  compared 
results  are  presented  in  Figure  42  and  Figure  43. 
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difference  between  data  and  McGinnis 


23  24  25  26  27  28  29  30 
Frequency,  GHz 


difference  between  data  and  McGinnis 


Figure  42.  (a)  Magnitude  Differences  of  Transmission  Coefficient  between  the 

Calculated  and  Simulated  Data  (b)  Phase  Differences  (Degrees)  of  Transmission 
Coefficient  between  the  Calculated  and  Simulated  Data 


difference  between  data  and  McGinnis  difference  between  data  and  McGinnis 


Frequency,  GHz  Frequency,  GHz 


Figure  43.  (a)  Magnitude  Differences  of  Reflection  Coefficient  between  the  Calculated 

and  Simulated  Data  (b)  Phase  Differences  (Degrees)  of  Reflection  Coefficient 
between  the  Calculated  and  Simulated  Data 
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According  to  the  Figures  42  and  43,  the  differences  between  the  calculated  and 
simulated  data  for  the  real  parts  of  the  transmission  and  reflection  coefficients  are  almost 
zero.  There  is  less  than  a  one  degree  phase  difference  in  both  the  transmission  and 
reflection  coefficients,  which  is  acceptable.  Therefore,  the  simulation  set  up  is  considered 
correct  and  the  results  are  acceptable. 

For  the  same  model  and  data  that  was  presented  above,  an  extraction  program  is 
used  to  retrieve  the  permittivity  and  permeability  from  the  scattering  parameters 
generated  from  MWS.  Figures  44  and  45  show  the  results  for  the  extracted  jur  and  sr 
using  MATLAB. 
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Figure  44.  Extracted  Real  Part  and  Imaginary  Part  of  the  Penneability  of  Model  5 
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Figure  45.  Extracted  Real  Part  and  Imaginary  Part  of  the  Permittivity  of  Model  5 


The  constitutive  parameters  were  successfully  extracted  from  the  simulated 
S-parameters.  Besides,  it  was  found  that  although  there  was  less  than  one  degree  phase 
difference  between  the  simulated  and  calculated  S-parameters,  the  differences  did  not 
have  a  big  influence  on  the  extracted  permeability  and  permittivity. 

2.  Model  6  -  “New”  Metamaterial — H  Type  Metamaterial  in  Waveguide 
In  order  to  compare  with  the  results  from  model  4,  which  used  the  free  space 
method  to  extract  the  constitutive  parameters,  the  dimensions  of  the  waveguide  were 
chosen  so  that  the  waveguide  can  operate  at  the  same  frequency  range  as  in  free  space 
(model  4).  The  related  dimensions  and  cut-off  frequencies  are  listed  in  Table  3  and  Table 
4. 
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U) 

m  =1 

m  =2 

m  =3 

n  =0 

9.37  GHz 

18.74  GHz 

28.11  GHz 

n  =1 

36.91  GHz 

40.32  GHz 

45.44  GHz 

n  =2 

72.01  GHz 

73.81  GHz 

76.73  GHz 

n  =3 

107.5  GHz 

108.72  GHz 

110.72  GHz 

Table  3.  Cut-off  Frequencies  for  Waveguide  with  a  =  16  mm  and  b  =  4.2  mm. 


(X) 

V  *'  c  /  mn 

m  =1 

m  =2 

m  =3 

n  =0 

6.25  GHz 

12.49  GHz 

18.74  GHz 

n  =1 

36.24  GHz 

37.82  GHz 

40.32  GHz 

n  =2 

71.67  GHz 

72.48  GHz 

73.81  GHz 

n  =3 

107.28  GHz 

107.82  GHz 

108.72  GHz 

Table  4.  Cut-off  Frequencies  for  Waveguide  with  a  =  24  mm  and  b  =  4.2  mm. 


According  to  Tables  2,  3,  and  4,  the  width  of  the  waveguide  a  can  only  be  chosen 
to  be  16  mm  so  that  the  simulation  can  run  in  the  frequency  range  of  14  to  18  GHz  and 
only  the  TEl0  mode  will  propagate  in  the  waveguide. 

The  schematic  of  the  H-type  inclusion  metamaterial  in  the  waveguide 
environment  is  shown  in  Figure  46.  The  boundary  condition  setting  is  the  same  as  in 
model  5.  The  waveguide  dimensions  are  a  =  16  nun  and  b  =  4.2  mm. 
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Figure  46.  Two  Unit  Cells  of  H  Type  Metamaterial  with  the  Same  Dimension  Described 

in  Model  4 

The  extracted  permeability  and  permittivity  are  shown  in  Figures  47  and  48.  Note 
that  there  is  a  second  spike  observed  in  both  the  real  and  imaginary  parts  of  permittivity. 
This  second  peak  only  occurs  in  s  when  the  material  is  in  the  waveguide  and  arises 
from  the  1/2 a  term  in  the  numerator  of  Eq.  (4. 1 8).  The  peak  disappears  as  a  — »  oo  and 
normal  incidence  is  simulated.  As  a  practical  matter,  the  material  would  not  be  used  near 
resonance  because  the  losses  are  high.  Most  published  papers  omit  data  near  the  resonant 
frequency,  a  range  known  to  be  hard  to  deal  with  in  the  retrieval  of  metamaterials 
parameters  [24]. 

One  possible  source  of  this  peak  is  that  the  extraction  formulas  assume  an 
isotropic  material.  Formulas  are  available  for  extracting  anisotropic  materials  in  free 
space,  but  not  for  waveguide  [23].  This  feature  needs  further  investigation. 
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Figure  47.  Extracted  Permeability  of  Model  6 


Figure  48.  Extracted  Permittivity  of  Model  6 


D.  SUMMARY 

In  this  chapter,  two  models  were  used  to  extract  the  penneability  and  permittivity 
from  the  simulated  S-parameters.  The  formulas  used  were  based  on  samples  in  a 
waveguide  environment,  and  therefore  the  image  principle  had  to  be  applied  in  order  to 
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treat  the  model  as  an  infinite  slab.  Test  cases  were  performed  to  verify  the  effectiveness  of 
the  simulation  set  up.  It  was  found  that  the  differences  between  the  simulated  and 
calculated  S-parameter  were  negligible.  Again,  the  resonant  frequency  of  the  metamaterial 
model  was  very  sensitive  to  the  mesh  cell  number;  therefore,  it  is  good  to  have  both 
simulation  and  measurement  data  to  compare.  The  models  have  shown  that  the  waveguide 
environment  program  can  work  not  only  for  nonnal  materials  but  also  metamaterials.  In 
comparing  model  4  with  model  6,  the  extracted  permeability  and  permittivity  were  very 
close.  This  result  showed  that  both  the  free  space  and  waveguide  method  can  be  used  to 
retrieve  the  constitutive  parameters  of  materials. 
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V.  CONCLUSION 


A.  SUMMARY  AND  CONCLUSIONS 

The  basic  concept  of  a  metamaterial  and  the  unusual  properties  of  DNG  materials 
were  discussed.  The  applications  of  metamaterials  in  waveguide  miniaturization, 
microwave  devices  and  other  possible  applications  were  presented  and  explained.  To 
employ  metamaterials  in  devices  it  is  necessary  to  know  their  effective  pennittivity  and 
permeability.  The  information  can  be  obtained  by  measurements  on  a  sample  in  free  space 
or  waveguide.  Alternately,  simulation  of  the  material  in  free  space  or  a  waveguide  could  be 
used. 

The  equations  used  in  the  free  space  enviromnent  to  extract  the  constitutive 
parameters  from  simulated  S-parameters  were  described  and  programmed  in  MATLAB. 
Test  cases  were  used  to  validate  the  correctness  of  the  simulation  settings.  Several  different 
types  of  materials,  including  the  normal  (homogeneous)  material  and  metamaterial  were 
simulated  in  MWS  and  a  retrieval  program  was  used  to  extract  the  pennittivity  and 
permeability  of  the  materials.  The  results  were  compared  with  published  data.  The  resonant 
frequency  of  the  results  and  its  relationship  to  the  dimensions  of  the  inclusions  and  number 
of  mesh  cells  at  the  simulation  starting  point  was  investigated  and  discussed. 

The  waveguide  enviromnent  retrieval  method  and  accompanying  equations  were 
discussed  in  Chapter  IV.  Boundary  conditions  were  properly  selected  so  that  a  small 
sample  in  the  waveguide  environment  could  be  treated  as  an  infinite  panel.  The  same 
metamaterial  (H-type  inclusions)  was  examined  using  both  the  free  space  and  waveguide 
methods  to  compare  the  extracted  results. 

It  was  shown  that  both  the  free  space  and  waveguide  methods  could  be  used  to 
successfully  retrieve  the  constitutive  parameters  from  the  scattering  parameters.  In  both 
free  space  and  waveguide  cases  the  number  of  mesh  cells  used  at  the  simulation  starting 
point,  the  dimensions  of  the  inclusions,  and  the  boundary  condition  settings  in  MWS  were 
the  key  parameters  for  successful  retrieval. 
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Finally,  the  simulation  and  extracted  results  of  the  effective  parameters  could  be 
used  as  data  for  designing  microwave  devices  with  metamaterials.  The  simulation  methods 
and  extraction  programs  make  the  effective  parameter  determination  much  more  efficient 
in  time  and  cost  than  building  large  panels  for  testing. 


B.  FUTURE  WORK 

In  this  thesis,  normal  and  DNG  materials  were  used  to  examine  the  retrieval 
programs.  Single  negative  materials  (SNG)  were  not  included.  Dealing  with  the  SNG 
materials  is  more  of  a  challenge  because  of  the  sign  ambiguity  in  the  square  root  of  //  •  s 
when  it  is  programmed.  More  inputs  must  be  acquired  to  decide  the  correct  sign  when 
using  SNG  materials.  Future  research  could  be  conducted  to  validate  the  programs  using 
thicker  materials  which  are  more  realistic  in  tenns  of  fabrication.  In  this  research  the 
sample  thickness  has  been  small  compared  to  wavelength.  This  condition  avoids  many 
ambiguities  that  occur  for  thick  samples.  Future  work  could  also  be  conducted  to  correct 
the  errors  due  to  the  nonzero  incidence  so  that  the  extracted  results  could  be  compared  with 
the  normal  incident  cases.  Finally,  the  cause  of  the  second  resonance  in  s  that  occurs  in 
the  waveguide  needs  to  be  examined  further. 
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APPENDIX 


A.  MATLAB  CODE  FOR  FREE  SPACE  ENVIRONMENT 

This  appendix  provides  the  MATLAB  code  used  to  retrieve  the  constitutive 
parameters  from  S-parameters  in  Chapter  III.  The  code  is  written  in  MATLAB  7.  This  code 
retrieves  and  plots  //  ,  sr ,  n,  and  Z  as  a  function  of  frequency  respectively. 


%###################################################### 

%  multiple  frequency  for  plane  wave  incident  on  a  panel 

%  sample  has  thickness  d,  free  space  case 

%###################################################### 

clear 

clc 

^~k  -k  ~k  ~k  -k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  -k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  -k  ~k  ~k  -k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  -k  ~k  ~k  ~k  -k  -k 

%  constants 
rad=pi/180; 
eps0=8 . 85e-12 ; 
mu0=4  *pi*le-7; 

etaO=sqrt (muO/epsO) ;  %  free  space 

c=l/ sqrt (epsO*muO) ; 

Q-'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

o 

%  variables 

d=2.5e-3;  %  thickness  of  the  sample 

fl=2e9;  df=18e6;  f2=20e9;  Nf=f loor ( (f 2-f 1 ) /df ) +1 ;  %  frequency  range  and 
spacing 

o 

%  load  raw  data 
load  Mlr.m; 
load  M2r.m; 
load  Plr.m; 
load  P2r.m; 

Magl=Mlr ( : , 2) ; 

Mag2=M2r ( : , 2 ) ; 

Phl=-unwrap (Plr ( : , 2) , 180) ; 

Ph2=-unwrap (P2r ( : , 2) , 180) ; 

9''k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

o 

%  freqyency  loop  begin 
for  it=l : Nf 

f=fl+ (it-1) *df ; 

Fr (it) =f ; 

Fghz (it) =f/le9; 
fghz=f /le9; 

Z0=eta0;  W=2*pi*Fr (it) ;  betaO=W/c; 

^•k  ~k  -k  -k  k:  -k  ~k  -k  ~k  -k  -k  -k  ~k  ~k  'k  'k  ~k  k:  ~k  -k  ~k  -k  ~k  ~k  -k  -k  ~k  -k  ~k  -k  ~k  -k  -k  -k  ^  ~k  -k  -k  ~k  -k  -k 

%  combine  magnitude  and  phase  of  the  Sll  and  S21 
R (it) =Magl (it) *exp ( j  *Phl (it) *rad)  ; 

T (it) =Mag2 (it) *exp ( j  *Ph2 (it) *rad)  ; 

9-'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

o 
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%  m  loop  begin 

Max=0;  %  - 

itt=0 ; 

for  m= -Max : Max 
itt=itt+l ; 


-m  value  set  up 


9-  *****  ; 


%  using  formulas  from  Robust  method 
Tp=T (it) ; 

Z21 (it) =sqrt ( ( (1+R(it) ) A2-TpA2) / ( (l-R(it) ) A2-TpA2) ) ; 

Z22 (it) =-sqrt ( ( (1+R(it) ) A2-TpA2) / ( (l-R(it) ) A2-TpA2) ) ; 
expinkdl (it) =Tp/ (l-R(it) * (Z21 (it) -1) / (Z21  (it) +1) ) ; 
expinkd2 (it) =Tp/ (l-R(it) * (Z22 (it) -1) / (Z22  (it) +1) ) ; 
if  abs (real (Z21 (it) ) ) >=0 . 005  &  real (Z21 (it) ) >=0 
expinkd (it) =expinkdl (it) ; 

Z2 (it) =Z2 1 (it) ; 

end 

if  abs (real (Z21 (it) ) ) >=0 . 005  &  real (Z21 (it) ) <0 
expinkd (it) =expinkd2 (it) ; 

Z2 (it) =Z22 (it) ; 

end 

if  abs (real (Z21 (it) ) ) <0 . 005  &  abs (expinkdl (it) ) <=1 
expinkd (it) =expinkdl (it) ; 

Z2 (it) =Z2 1 (it) ; 

end 

if  abs (real (Z21 (it) ) ) <0 . 005  &  abs (expinkdl (it) ) >1 
expinkd (it) =expinkd2 (it) ; 

Z2 (it) =Z22 (it) ; 

end 

ni (it)  =-1/ (beta0*d) *j *real (log (expinkd (it) )) ;  %  imag  of  n 

nr (it,  itt) =1/ (beta0*d) * (imag (log (expinkd (it) )) +2*m*pi) ;  %real  of  n 

n2 (it,  itt) =nr (it, itt) +ni  (it) ; 

Er2 (it, itt) =n2 (it, itt) / Z2  (it)  ; 

Mr 2 (it, itt) =n2 (it, itt) *Z2 (it) ; 
end  %  end  of  m  loop 
end  %  end  of  frequency  loop 

figure ( 1 ) 

subplot  (121)  ; 

plot (Fghz,  real (Mr2) ) ; 

xlabel (' Frequency  in  GHz') 

ylabel ( ' Re ( \mu) ' ) 

subplot  ( 122 ) 

plot (Fghz, imag (Mr2) ) ; 

xlabel (' Frequency  in  GHz') 

ylabel ( ' Im ( \mu) ' ) 


figure (2 ) 

subplot (121) ; 

plot (Fghz, real (Er2) ) ; 

xlabel ( ' Frequency  in  GHz ' 

ylabel (' Re ( \epsilon) ') 

subplot  ( 122 ) 

plot (Fghz, imag (Er2) ) ; 

xlabel ( ' Frequency  in  GHz ' 

ylabel (' Im ( \epsilon) ') 
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figure ( 3 ) 

subplot (121) ; 

plot (Fghz, real (n2) ) ; 

xlabel ( ' Frequency  in  GHz ' ) 

ylabel ( ' Re (n) ' ) 

subplot  ( 122 ) 

plot (Fghz, imag (n2) ) ; 

xlabel ( ' Frequency  in  GHz ' ) 

ylabel ( ' Im (n) ' ) 

figure ( 4 ) 

subplot  (121)  ; 

plot (Fghz, real  (Z2) )  ; 

xlabel ( ' Frequency  in  GHz ' ) 

ylabel ( ' Re ( Z ) ' ) 

subplot  ( 122 ) 

plot (Fghz, imag ( Z  2 )  ) ; 

xlabel ( ' Frequency  in  GHz ' ) 

ylabel ( ' Im ( Z ) ' ) 
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B.  MATLAB  CODE  FOR  WAVEGUIDE  ENVIRONMENT 

This  appendix  provides  the  MATLAB  code  used  to  retrieve  the  constitutive 
parameters  from  S-parameters  in  Chapter  IV.  The  code  is  written  in  MATLAB  7.  This 
code  retrieves  and  plots  //,. ,  sr ,  n,  and  Z  as  a  function  of  frequency  respectively. 


%####################################################################### 

%  multiple  frequency  in  a  waveguide  environment 

%  sample  has  thickness  delta,  total  length  of  the  waveguide  is  dl+delta+d2 

%####################################################################### 

clear 

clc 

a=l 6e-3 ; 

delta=2 .  274e-3; 

dl=le-3 ; 

d2=le-3 ; 

eps0=8 . 85e-12  ; 

mu0=4*pi*le-7; 

c=l/sqrt (epsO*muO) ; 

rad=pi/180; 

Z 0=3 7 7 ; 

Q-'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

o 

load  Mlr.m; 
load  M2r.m; 
load  Plr.m; 
load  P2r.m; 

Magl=Mlr ( : , 2)  ; 

Mag2=M2r ( : , 2 ) ; 

Phl=-unwrap (Plr ( : , 2) , 180) ; 

Ph2=-unwrap ( P2r ( : , 2 )  ,180)  ; 

★  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  ~k  -k  -k  -k  ★  -k  -k  -k  -k  -k  -k  -k  ~k  -k  -k  -k  ~k  -k  -k  -k  ★ 

it=0  ; 

f l=14e9;  df=4e6;  f2=18e9;  Nf=floor ( (f2-fl) /df ) +1; 
for  it=l : Nf 

f=fl+ ( it— 1 ) *df ;  %  must  higher  than  the  cut-off  frequency  in  free  space 

Fr (it) =f ; 

w=2*pi*Fr (it) ;  wave=c/f; 

Rm (it) =Magl (it) *exp ( j  *Phl (it) *rad) ; 

Tm (it) =Mag2 (it) *exp ( j  *Ph2 (it) *rad) ; 
betaO (it) =sqrt ( (w/c) A2- (pi/a) A2) ; 

end 

9''k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

o 

%  make  use  of  equation  from  Thomas 

nf req=it; 

for  it=l:nfreq 

Tt (it) =Tm (it) *exp (-j*beta0 (it) * (dl+d2 ) ) ;  % - T 

Rr ( it ) =Rm (it) *exp (- j  *2*beta0  (it ) *dl ) ;  % - R 

Q-'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

%  to  get  Z2 

b= (Rr (it) A2-Tt (it) A2+l) / (2*Rr (it) ) ; 

R21=b+ (sqrt (bA2-l) ) ; 

R22=b- (sqrt (bA2-l) ) ; 
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Rn=nan; 

if  abs (R2 1 ) <=1 

RR=R2 1 ;  Rn=RR;  % - R2 

end 

if  abs (R22 ) <=1 

RR=R2  2 ;  Rn=RR; 

end 

if  (abs(R21)<l  &  abs(R22)<l)  |  (abs(R21)==l  &  abs(R22)==l) 

disp(['both  solutions  give  |R|<1  for  fghz  =  ' , num2str (fghz) ] ) 

end 

Z= (1+RR) / (1-RR)  ; 

sgn=l;  if  real(Z)<0,  sgn=-l;  end 
Z2 (it) =Z*sgn; 

o, _ 

o 

k=0;  %  need  to  be  determined 

o,  _ 

o 

K=Rr (it) A2-Tt (it) A2+l; 
rl= (K+sqrt (KA2-4*Rr (it) A2) ) /2/Rr (it)  ; 
r2= (K-sqrt (KA2-4*Rr (it) A2) ) /2/Rr (it) ; 
if  abs (rl ) <=1 
r=rl  ; 

end 

if  abs (r2 ) <=1 
r=r2  ; 

end 

t=  (Rr  (it)  +Tt  (it)  -r)  /  (1-r*  (Rr  (it)  +Tt  (it)  )  )  ; 
gamma= ( j  *2*k*pi-log (t) ) /delta; 

gammaO=- j  *2 *pi*sqrt (Fr (it) A2 *muO*epsO- (1/2/a) A2) ; 

Mr (it) =muO* (gamma* (1+r) /gammaO/ (1-r) ) ; 

Er (it) = ( ( (1/2/a) A2- (gamma/ 2 /pi) A2) /Mr (it) /Fr (it) A2) ; 

Mre (it) =Mr (it) /muO; 

Ere (it) =Er (it) /epsO; 
n2 (it) =Mre (it) / Z2  (it)  ; 
end  %  end  of  frequency  loop 

figure ( 1 ) 
subplot  (121)  ; 
plot (Fr/le9, real (Mre) ) ; 
xlabel (' Frequency  in  GHz') 
ylabel ( ' Re ( \mu) ' ) 
subplot  ( 122 ) 
plot (Fr/le9, imag (Mre) ) ; 
xlabel (' Frequency  in  GHz') 
ylabel ( ' Im ( \mu) ' ) 

figure  (2 ) 
subplot  (121)  ; 
plot (Fr/le9, real (Ere) ) ; 
xlabel (' Frequency  in  GHz') 
ylabel (' Re ( \epsilon) ') 
subplot  ( 122 ) 
plot (Fr/le9, imag (Ere) ) ; 
xlabel (' Frequency  in  GHz') 
ylabel (' Im ( \epsilon) ') 
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figure ( 3 ) 

subplot (121) ; 

plot (Fr/le9, real (n2) ) ; 

xlabel ( ' Frequency  in  GHz ' ) 

ylabel ( ' Re (n) ' ) 

subplot  ( 122 ) 

plot (Fr/le9, imag (n2) ) ; 

xlabel ( ' Frequency  in  GHz ' ) 

ylabel ( ' Im(n) ' ) 

figure ( 4 ) 

subplot  (121)  ; 

plot (Fr/le9, real  (Z2) )  ; 

xlabel ( ' Frequency  in  GHz ' ) 

ylabel ( ' Re ( Z ) ' ) 

subplot  ( 122 ) 

plot (Fr/le9, imag (Z2) ) ; 

xlabel ( 1  Frequency  in  GHz ' ) 

ylabel ( ' Im ( Z ) ' ) 
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